The macroevolutionary singularity of snakes.

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Repeated evolution of blanched coloration in a lizard across independent white-sand habitats.

The White Sands lizards of New Mexico are a rare and classic example of convergent evolution where three species have evolved blanched coloration on the white gypsum dunes. Until now, no geological replicate of the pattern had been described. However, one of the White Sands species, the lesser earless lizard (Holbrookia maculata), has been discovered to also inhabit the Salt Basin Dunes of Texas, where it has also evolved a blanched morph. We here present a first phenotypic and genetic description of the Salt Basin Dunes population of H. maculata. Phylogenetic inference based on a housekeeping gene (ND4) and a classic candidate gene in the melanin-synthesis pathway (Melanocortin 1 Receptor; Mc1r) shows the newly discovered population as an independent lineage, with no evidence of genetic parallelism in the coding region of Mc1r. Initial morphological data suggest that while this population displays convergent evolution in blanched coloration, there are divergent patterns in limb length and habitat use behavior between the gypsum environments. Our findings present the White Sands/Salt Basin Dunes as an exceptionally promising comparative model for studies of adaptation and convergent evolution.

Genetic and Ecogeographic Controls on Species Cohesion in Australia's Most Diverse Lizard Radiation.

AbstractSpecies vary extensively in geographic range size and climatic niche breadth. If range limits are primarily determined by climatic factors, species with broad climatic tolerances and those that track geographically widespread climates should have large ranges. However, large ranges might increase the probability of population fragmentation and adaptive divergence, potentially decoupling climatic niche breadth and range size. Conversely, ecological generalism in large-ranged species might lead to higher gene flow across climatic transitions, increasing species' cohesion and thus decreasing genetic isolation by distance (IBD). Focusing on Australia's iconic Ctenotus lizard radiation, we ask whether species range size scales with climatic niche breadth and the degree of population isolation. To this end, we infer independently evolving operational taxonomic units (OTUs), their geographic and climatic ranges, and the strength of IBD within OTUs based on genome-wide loci from 722 individuals spanning 75 taxa. Large-ranged OTUs were common and had broader climatic niches than small-ranged OTUs; thus, large ranges do not appear to simply result from passive tracking of widespread climatic zones. OTUs with larger ranges and broader climatic niches showed relatively weaker IBD, suggesting that large-ranged species might possess intrinsic attributes that facilitate genetic cohesion across large distances and varied climates. By influencing population divergence and persistence, traits that affect species cohesion may play a central role in large-scale patterns of diversification and species richness.

No link between population isolation and speciation rate in squamate reptiles.

Rates of species formation vary widely across the tree of life and contribute to massive disparities in species richness among clades. This variation can emerge from differences in metapopulation-level processes that affect the rates at which lineages diverge, persist, and evolve reproductive barriers and ecological differentiation. For example, populations that evolve reproductive barriers quickly should form new species at faster rates than populations that acquire reproductive barriers more slowly. This expectation implicitly links microevolutionary processes (the evolution of populations) and macroevolutionary patterns (the profound disparity in speciation rate across taxa). Here, leveraging extensive field sampling from the Neotropical Cerrado biome in a biogeographically controlled natural experiment, we test the role of an important microevolutionary process-the propensity for population isolation-as a control on speciation rate in lizards and snakes. By quantifying population genomic structure across a set of codistributed taxa with extensive and phylogenetically independent variation in speciation rate, we show that broad-scale patterns of species formation are decoupled from demographic and genetic processes that promote the formation of population isolates. Population isolation is likely a critical stage of speciation for many taxa, but our results suggest that interspecific variability in the propensity for isolation has little influence on speciation rates. These results suggest that other stages of speciation-including the rate at which reproductive barriers evolve and the extent to which newly formed populations persist-are likely to play a larger role than population isolation in controlling speciation rate variation in squamates.

Congruence and Conflict in the Higher-Level Phylogenetics of Squamate Reptiles: An Expanded Phylogenomic Perspective.

Genome-scale data have the potential to clarify phylogenetic relationships across the tree of life but have also revealed extensive gene tree conflict. This seeming paradox, whereby larger data sets both increase statistical confidence and uncover significant discordance, suggests that understanding sources of conflict is important for accurate reconstruction of evolutionary history. We explore this paradox in squamate reptiles, the vertebrate clade comprising lizards, snakes, and amphisbaenians. We collected an average of 5103 loci for 91 species of squamates that span higher-level diversity within the clade, which we augmented with publicly available sequences for an additional 17 taxa. Using a locus-by-locus approach, we evaluated support for alternative topologies at 17 contentious nodes in the phylogeny. We identified shared properties of conflicting loci, finding that rate and compositional heterogeneity drives discordance between gene trees and species tree and that conflicting loci rarely overlap across contentious nodes. Finally, by comparing our tests of nodal conflict to previous phylogenomic studies, we confidently resolve 9 of the 17 problematic nodes. We suggest this locus-by-locus and node-by-node approach can build consensus on which topological resolutions remain uncertain in phylogenomic studies of other contentious groups. [Anchored hybrid enrichment (AHE); gene tree conflict; molecular evolution; phylogenomic concordance; target capture; ultraconserved elements (UCE).].

Is genomic diversity a useful proxy for census population size? Evidence from a species-rich community of desert lizards.

Species abundance data are critical for testing ecological theory, but obtaining accurate empirical estimates for many taxa is challenging. Proxies for species abundance can help researchers circumvent time and cost constraints that are prohibitive for long-term sampling. Under simple demographic models, genetic diversity is expected to correlate with census size, such that genome-wide heterozygosity may provide a surrogate measure of species abundance. We tested whether nucleotide diversity is correlated with long-term estimates of abundance, occupancy and degree of ecological specialization in a diverse lizard community from arid Australia. Using targeted sequence capture, we obtained estimates of genomic diversity from 30 species of lizards, recovering an average of 5,066 loci covering 3.6 Mb of DNA sequence per individual. We compared measures of individual heterozygosity to a metric of habitat specialization to investigate whether ecological preference exerts a measurable effect on genetic diversity. We find that heterozygosity is significantly correlated with species abundance and occupancy, but not habitat specialization. Demonstrating the power of genomic sampling, the correlation between heterozygosity and abundance/occupancy emerged from considering just one or two individuals per species. However, genetic diversity does no better at predicting abundance than a single day of traditional sampling in this community. We conclude that genetic diversity is a useful proxy for regional-scale species abundance and occupancy, but a large amount of unexplained variation in heterozygosity suggests additional constraints or a failure of ecological sampling to adequately capture variation in true population size.

Does Population Structure Predict the Rate of Speciation? A Comparative Test across Australia's Most Diverse Vertebrate Radiation.

Population divergence is the first step in allopatric speciation, as has long been recognized in both theoretical models of speciation and empirical explorations of natural systems. All else being equal, lineages with substantial population differentiation should form new species more quickly than lineages that maintain range-wide genetic cohesion through high levels of gene flow. However, there have been few direct tests of the extent to which population differentiation predicts speciation rates as measured on phylogenetic trees. Here, we explicitly test the links between organismal traits, population-level processes, and phylogenetic speciation rates across a diverse clade of Australian lizards that shows remarkable variation in speciation rate. Using genome-wide double digest restriction site-associated DNA data from 892 individuals, we generated a comparative data set on isolation by distance and population differentiation across 104 putative species-level lineages (operational taxonomic units). We find that species show substantial variation in the extent of population differentiation, and this variation is predicted by organismal traits that are thought to be proxies for dispersal and deme size. However, variation in population structure does not predict variation in speciation rate. Our results suggest that population differentiation is not the rate-limiting step in species formation and that other ecological and historical factors are primary determinants of speciation rates at macroevolutionary scales.

Predator Perspective Drives Geographic Variation in Frequency-Dependent Polymorphism.

Color polymorphism in natural populations can manifest as a striking patchwork of phenotypes in space, with neighboring populations characterized by dramatic differences in morph composition. These geographic mosaics can be challenging to explain in the absence of localized selection because they are unlikely to result from simple isolation-by-distance or clinal variation in selective regimes. To identify processes that can lead to the formation of geographic mosaics, we developed a simulation-based model to explore the influence of predator perspective, selection, migration, and genetic linkage of color loci on allele frequencies in polymorphic populations over space and time. Using simulated populations inspired by the biology of Heliconius longwing butterflies, Cepaea land snails, Oophaga poison frogs, and Sonora ground snakes, we found that the relative sizes of predator and prey home ranges can produce large differences in morph composition between neighboring populations under both positive and negative frequency-dependent selection. We also demonstrated the importance of the interaction of predator perspective with the type of frequency dependence and localized directional selection across migration and selection intensities. Our results show that regional-scale predation can promote the formation of phenotypic mosaics in prey species, without the need to invoke spatial variation in selective regimes. We suggest that predator behavior can play an important and underappreciated role in the formation and maintenance of geographic mosaics in polymorphic species.

Stable isotope ecology of a hyper-diverse community of scincid lizards from arid Australia.

We assessed the utility of stable isotope analysis as a tool for understanding community ecological structure in a species-rich clade of scincid lizards from one of the world's most diverse lizard communities. Using a phylogenetic comparative framework, we tested whether δ15N and δ13C isotopic composition from individual lizards was correlated with species-specific estimates of diet and habitat use. We find that species are highly divergent in isotopic composition with significant correlations to habitat use, but this relationship shows no phylogenetic signal. Isotopic composition corresponds to empirical observations of diet for some species but much variation remains unexplained. We demonstrate the importance of using a multianalytical approach to questions of long-term dietary preference, and suggest that the use of stable isotopes in combination with stomach content analysis and empirical data on habitat use can potentially reveal patterns in ecological traits at finer scales with important implications for community structuring.